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PREFACE TO PART IV. 

The first part of this work which contains prob- 
lems in machine design was pubUshed in 1912, and the 
favorable reception accorded it gave the author the 
courage to publish the rest of the problems which he 
had collected and used in his classes for some years. 
The second and third parts contain problems selected 
from various branches of hydraulics and thermody- 
namics respectively, in the solution of which it is neces- 
sary to use calculus and analytic geometry. The part 
now offered contains problems selected from various 
topics in mechanics of materials, and the last part con- 
tains problems on electrical engineering. A student or 
an engineer who wishes to review calculus or analytics, 
or to acquire faciUty in appUcations of higher mathe- 
matics to engineering problems, may thus select at first 
the part of the work which deals with problems with 
which he is most f amiUar, or in which he is particularly 
interested. 

The reader will find the author's views on teaching 
mathematics to engineering students in the Preface to 
Part I and in the Dialogue following that preface. He 
is also referred to Part I for a Ust of reference works on 
mathematics and for an Appendix entitled "What a 
Senior in Engineering ought to know about Mathe- 
matics." 

The author wishes to acknowledge gratefully the 
assistance of Mr. A. C. Stevens, M.E., mstructor m 
Cornell University, who read the proofs and made 
some valuable suggestions for the text. 

Cornell University, Ithaca, N. Y. 

ui 
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PART IV. 

PROBLEMS ON 
MECHANICS OF MATERIALS. 



ENGINEERING MATHEMATICS. 



CHAPTER I. 

STRENGTH OF SIMPLE HORIZONTAL BEAMS. 

Let a straight horizontal beam of uniform cross 
section (Fig. 1) be freely supported on two knife- 
edges, and loaded with vertical forces or weights. 




Fig. 1. — A loaded horizontal beam. 

Some of these forces or loads may be concentrated, 
such as P and Q, others uniformly or non-uniformly 
distributed along the whole span or over a portion of it. 
If a load is distributed it may be given as q kg. per 
running meter of the span, where q may be either a 
constant quantity or a function of the distance x from 
one of the supports. 

Under the influence of these forces the beam is de- 
flected, and its longitudinal axis AB becomes a curve. 

3 
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The cross sections of the beam which are originaUy per- 
pendicular to this axis, and which he in paraUel planes 
before .the forces are appUed, are now deformed by 
the stresses. A fundamental assumption is made in the 
theory of the strength of materials that as long as the 






Fig. 2. — The axis and two plane cross sections of a loaded beam. 

stresses in the beam do not exceed a certain " elastic ." 
limit, the plane cross sections of the beam remain plane 
in their new distorted position (Fig. 2). The beam 
may be considered to consist of longitudinal fibers 
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parallel to the axis, and the deformation of the beam 
causes a longitudinal compression along the fibers on 
the concave side of the beam (66') and a longitudinal 
tension along the fibers on the convex side (cc') . There 
is a set of fibers, such as aa^, which do not change in 
length during the deformation. They all lie before 
the deformation in a horizontal plane called the neutral 
plane. All the fibers above this plane are in a state 
of compression, those below the neutral plane are in 
a state of tension. 

It wUl be seen from Fig. 2 that maximum longitudmal 
stresses take place in the extreme fibers 66' and cc^, 
so that if the stresses in these fibers do not exceed a 
safe permissible hmit, the whole beam is safe.* The 
maximum stress varies from cross section to cross sec- 
tion of the beam, depending upon the position of the 
appUed forces. Since the action of these forces is to 
bend the beam, the internal stresses in a cross section 
are proportional to the bending moment of the appUed 
forces with respect to that cross section. For instance, 
in Fig. 1, the bending moment M with respect to cross 

section N is 

M = Bx- Qv, 

where B is the force reaction of the right-hand support. 
In this equation force Q is taken with the minus sign 
because its moment is opposite to that of force B. 
Since the beam is in equilibrium, the numerical value 
of the moment of the forces to the left of the cross 
section N is also equal to Af , for otherwise the beam 
would be turning imder the influence of an imbalanced 

* The beam must also be checked for the shearing stress near the 
neutral plane. This stress is not considered here, and is assumed to 
be small; otherwise, the maximum combined stress does not necessarily 
take place in the extreme fibers. 



6 ENGINEERING MATHEMATICS. [Chap. I. 

moment. We took the moment of the forces to the 
right of N because in this case its expression happens 
to be simpler. 

Let the distance from the neutral axis (Fig. 2) to 
the extreme compressed fiber be zi cm., and let the 
stress in this fiber be pi kg, per sq. cm. The relation- 
ship between this stress and the bending moment M is 
M = pi//2i (1) 



Pio. 3. — The cross eection of a beam. 

In this formula / stands for the moment of inertia of 
the cross section with respect to the horizontal axis 
which hes in its plane and passes through the center 
of gravity of the cross section.* Thus, referring to 
Fig. 3, which represents a cross section of the beam 

.... (2) 



'^r/ 



* The student will find proof of this fundamental formula ii 
standard work on mechanics of materials. 
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where ds is an infinitesimal horizontal strip of the area, 
and z is the distance of this strip from the horizontal 
axis passing through the center of gravity. 

For the Tna.ximiim tensile stress in a cross section we 
can write, by analogy with Eq. (1), and using the nota- 
tion in Fig. 2 : 

M = ThI/zi . (3) 

Prob. I. — Check the physical dimensions in formula (1). 

SohUion, — 

ke 

kg. X cm. = — ^ X cm.*/cm., 

cm.2 

or, after simplification, 

kg. X cm. = kg. X cm., 
which is an identity. 

Prob. 2. — Prove that for a rectangular cross section of width h 
and height A the moment of inertia/ = ^bh^,]a,ndl/zi = I/Z2 = ihh^. 

Prob. 3. — Referring to the preceding problem, a rectangular 
beam is to be sawed out of a log of diameter D. Determine 6 and 
h in such a way as to get a beam that will carry the greatest 
possible load. Hint. J hh^ = max., and 6^ + h^ = D* = const. 

Prob. 4. — A log is to be sawed into thin boards of equal width, 
and it is desired to utilize the material so as to have as Uttle waste 
as possible. Compare per cent weight in the log utilized for boards 
with that which could be utiUzed for a beam of maximum strength 
as required in the preceding problem. 

Prob. 5. — With a given factor of safety, how much more load 
can the original log carry than the strongest rectangular beam cut 
from it? 

PiQb. 6. — It is proved in mechanics of materials that the deflec- 
tion of a beam under a given load is proportional to the moment of 
inertia /, and not to the ratio I/zi, Modify Prob. 3 to read that 
the beam must be the stiffest possible and not the strongest possible. 

Prob. 7. — Make a sketch showing the cross section of a log 
and of three beams that may be sawed out of it, namely, the beam 
of mftYiTniiTn cross scction, the strongest beam, and the stiffest 
beam. 

Prob. 8. — Make a sketch of the cross section of an I beam, 
denote its dimensions by suitable letters and deduce a general 
expression for its moment of inertia. 
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Piob. 9. — A rectangular horizontal plate (Fig. 4) is fastened in 
a wall and is loaded with a vertical force P at one edge. What is 
the weakest plane along which the fracture is most liable to occur? 
Solution. — Let tbe wavy line represent a possible plane of failure 
determined by an unknown angle a. Applyii^ Eq. (I) we get 
(see Prob. 2) Px = p ■iyh\ or p = const. X (x/y). Thus, the 
weakest plane, viz., the one for which the stress p is a maximum, 
is determined by the condition x/y = max. Both x and y are 
functions of angle a, and we have x = cama,y =^ c/cos a, so that 

sin a cos a = maX. But 2 sin « cos a = sin 2 a, and its muTirrninn 

value is reached when a = 45°.' 



Fia. 4. — A horizontal plate loaded at one comer. 

Prob. 10. — Find by direct differentiation the mamnum value 
of sin a cos « required in the preceding problem. 

Prob. II. — Referring to Prob. 9, how much stronger is the 
plate along planes corresponding to a = 40° and « = 50° than it is 
for the 45-degree plane? 

Prob. II. — Show that the answer to Prob. 9 is correct only 
when l>2c (Fig. 4). When i = 2 c the plate is theoretically 

* The conditions in this problem are approximately those which 
obtain between two meshing gears in the moat unfavorable case when 
the pressure is transmitted at one comer of the tooth. The tooth ia 
most Uable to break abng a 45-d(^ree plane. 



Chap. I.] STRENGTH OF SIMPLE HORIZONTAL BEAMS, 9 

just as liable to break diagonally along a = 45° as along the plane 
AB. When I <2c the plane AB is the weakest cross section. 

Prob. 13. — Prove that the plate (Fig. 4) is strongest along the 
diagonal plane AC. 

Prob. 14. — Find a proof of formula (1) in some book, and give 
it in your own words. 

Let the beam {Fig. 1) be loaded with the weight Q 
only, at a distance u from the right-hand support, and 
let it be required to find the dangerous cross section of 
the beam, that is, the one for which the bending moment 
due to the force Q is a maximum. Let this cross sec- 
tion be Ny at an imknown distance x from the support 
B. According to the rule for the addition of parallel 
forces, the reactions of the supports are inversely as 
the distances from the load Q, and we have 

Reaction A = jQ, (4) 

Reaction B = — = — Q, . . . (5) 

the sum of the two reactions being equal and opposite 
to the force Q. 
The. bending moment with respect to point N is 

M^^Q(l-x), (6) 

and this expression holds true as long as no force is 
appUed between A and N* The moment M reaches 
its maximum for the smallest possible value of x, or. 
in this case when x = u. Thus, with one concentrated 

* In this chapter x generally denotes the distance from B to some 
point on the beam, such as iV; u denotes the distance from B to the 
point of application of a force. In different problems either x or, i» 
or both may be variable. 



J J 
J J 
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load the maximuin value of the bending moment takes 
place mider the load, and is equal to 

M = |(Z-x)Q (7) 

Eq. (6) shows that with a fixed load {u = const.) the 
bending moment varies along the beam according to a 
straight-Une law. This is shown m Fig. 1 by the 
broken line ACB. 

Prob. 15. — Show that with a moving load Q, the bending 
moment is a maximum when the load is in the center of the span. 

Prob. 16. — A uniformly distributed load of q kg. per meter of 
span is appUed to the beam as shown in Fig. 1, viz., between the 
abscissae xz and Xi. Calculate the bending moment due to this 
load at a given point A"^. Solvlion. — The infinitesimal load at a 
distance u from B\&q dUy so that the moment at N due to this load 
is, according to Eq. (6), dM = u(J — x) q du/L Integrating this 
expression between the limits u == Xz and u = 0:4, the desired moment 
due to the whole load is obtained. 

Prob. 17. — The beam is loaded simultaneously with a concen- 
trated weight Q at a distance u from support B (Fig. 1) and with 
a distributed load q per unit length between xz and X4. Find the 
abscissa x for which the total bending moment is a maximum. 

Prob. 18. — A load is non-uniformly distributed over a portion 
of the span, so that q = / (w). Write down the general expression 
for M at point N, and apply it to a few simple cases, for example, (a) 
the load increases or decreases according to a straight-line law from 
the supports to the center of the span; (6) the load varies according 
to a parabohc law, either over the whole span, or from both supports 
towards the center; (c) the load varies according to an exponential 
law; (d) the load varies inversely as the nth power of the distance 
from one of the supports; (e) the load varies according to the 
sine law. In each case find point N for which the bending moment 
is a maximum. 

Prob. 19. — Two unequal concentrated loads, P and Q (Fig. 1), 

are rigidly connected together at a constant distance a, and move 

^long the span. Find the greatest bending moment which this 

t iyi^tem of loads exerts upon the span, and the dangerous cross sec- 
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tion of the beam. Two connected moving loads may be due in 
practice to two pairs of wheels of a truck moving along a bridge. 
Solution. — The distribution of moment due to Q is represented by 
the line ACB, that due to P by the Une ADB. The combined mo- 
ment is obtained by adding the ordinates of these two curves. It 
will thus be seen that the maximum moment is under one of the 
loads, and not between or outside them. Assume first that max. M 
is under Q. Then for any position of Q, determined by an abscissa 
u = Xy the total moment is M = Qx {I — x) /I + P {x + a){l — x)/ly. 
according to Eqs. (7) and (6). . Take a derivative of M with re- 
spect to X and equate it to zero. The value of x from this equa- 
tion will give the position of the dangerous cross section. Make 
similar calculations for weight P, and compare the values of max. M 
so obtained. 

Prob. 20. — Referring to the preceding problem, prove that with 
P equal to Q, the combined load exerts its maximum moment when 
the distances from the center of the span to P and Q are equal to i a 
and 1 a respectively. 

Prob. 21. — Extend Prob. 19 to the case of three unequal weights, 
P, Q, and /?, at distances a and h respectively. 

Prob. 22. — Extend Prob. 19 to the case when, in addition to P 
and Qj the beam is subjected to a uniform load q, for instance, its 
own weight. Note that in this case the maximum moment is not 
necessarily under one of the concentrated loads. 



CHAPTER II. 

THE ELASTIC CURVE OF A BEAM. 

An elastic curve is the curve which represents the 
shape of the neutral axis of a beam (Fig. 1) after flexure. 
It is necessary in some practical problems to know the 
equation of this curve ; for instance, 

(a) To calculate the maximiun deflection of the 

beam; 
(6) To determine the strength of a " constrained " 

beam, or a beam fixed at the ends, 
(c) In the theory of continuous girders, or beams 

resting on more than two supports. 

Let the abscissse x of the elastic curve be measured 
from the right-hand support, and let the ordinates y 
be considered positive downwards. It is proved in 
the mechanics of materials that the differential equa- 
tion of the elastic curve is 

dx^ EI' ^ ^ 

where M and / have the same meaning as in the pre- 
ceding chapter, and E is the modulus of elasticity or a 
coefficient which characterizes the extensibifity of the 
material. If the load on the beam is given, M is a 
known function of x, so that Eq. (1) may be integrated 
twice and the ordinate y of the elastic curve expressed 
as a fxmction of x. 

The minus sign is necessary in Eq. (1) if the moment 
M is to be considered positive when it bends the beam 

12 
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in such a way that the upper side is concave. An 
inspection of Fig. 1 will show that in this case dhf/dx^ is 
negative for the whole curve, and in order to make it 
positive the minus sign is necessary. 

A few examples will make the use of Eq. (1) clear. 

Case A. A beam is loaded with a concentrated 
weight P at the center of the span. The reaction of 
each support is | P, and the bending moment at a point 
determined by the abscissa x, when x <\l,'^M =\ Px. 
Substituting this value in Eq. (1) we get 

d^ _ P 
dx^ 2EI^' 

Integrating this expression once, we get 

d^^2^^^'"*^'^' .... (2) 

where Ci is the constant of integration. Integrating, 
again, gives 

y^^j{C, + Crx--\7?). ... (3) 

The constants Ci and C2 are determined from the follow- 
ing conditions : when 

"" ^^' dx "' 
when X =0, y = 0. 

Using the first condition in Eq. (2), we get 

the second condition used in Eq. (3) gives 

Co = 0, 
so that finally 

y=^(lPx-l:i?) (4) 
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The deflection is a maximum in the center of the span, 
and its value is 

^'""■"482?/ ^^ 

Case B. Let now the weight P be appUed at a 
distance c from the right-hand support. In this case, 
on account of the imsymmetrical load, it is necessary to 
consider parts of the beam to the left and to the right 
of the load, because it is not possible to say offhand 
on which side the deflection is a maximum. The re- 
action of the support B is equal to P {I — c) /l, so that 
by analogy with Eq. (6) in the preceding chapter, 
the moment at a distance x from B (for x < c) is 
M = Px {I — c) /I Thus for a point on the right-hand 
side of the beam we have : 

__ d^ _ Pjl-c) ,^. 

dx^~ Ell ""' ^^^ 

or, after integration. 



and 



I = ^1-^ (C - i^). (7) 



y-^^^J^^(C'+Cz-^^),. . . (8) 



where C and C are constants of integration, to be 
determined from the terminal conditions of the prob- 
lem; namely, when x = 0, t/ = 0, so that Eq. (8) 
gives C = 0. The constant C must be determined 
from the condition that the deflection y is the same for 
the right-hand side of the beam as for the left-hand 
side, when x = c. Thus, it is necessary to write an 
equation corresponding to Eq. (8), for the left-hand 
side of the beam. 
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Were we to use the same origin B for the left-hand 
side of the beam this equation would be of a different 
form from Eq. (8), and the determination of the con- 
stants of integration would be somewhat tedious. In 
order to remedy this, we transfer the origm to A, and 
write directly, by analogy with Eqs. (7) and (8) : 

i-^^(''-^"' <»> 

and 

y=^^^^{Tt-\^). . . . (10) 

In these equations t is the abscissa measured from the 

origm A, r is a constant of mtegration, and Ci is the 

distance from the support A to the force P (Fig. 1), 

so that 

c + Ci = Z (11) 

At X = c, and t — Ci, the two elastic curves, repre- 
sented by Eqs. (8) and (10), have a common ordinate 
and a common tangent. Or 



\dxh u L: 



(11a) 



Using these conditions and Eq. (11), we get 

ci {Cc-k(^)=c (rci - i ci»), . . . (12) 
- ci (C - ^ c») = c (r - I ci«). . . . (13) 

Solving these as simultaneous equations for C and F, 

we obtain 

C = ic(c+2ci), 

r = ici(ci+2c).j 

Substituting these values in Eqs. (8) and (10), the 
equations of the two parts of the elastic curve are ob- 
tained in their final form. 



(14) 
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To find the maximum deflection, it is necessary to 
equate dy/dx and dy/dt to zero, because the deflection 
may reach its maximum on either side of the applied 
load. Equating expression (7) to zero we get for the 
abscissa Xm at which the deflection is a maximmn: 



x. = V2C = v/^^^^. . . (15) 
Similarly, equating expression (9) to zero gives 



t^ = V2T = ^ '^^'^+^'\ . . 



(16) 



Let c < Ci, that is, let the load be applied at the right 
of the center of the span. It can be shown that in 
this case Xm > c, and tm < Ci, consequently, * Eq. (15) 
has no physical meaning, and the actual maximum 
takes place at the abscissa tm determined by Eq. (16). 
On the other hand, if c > Ci, the deflection is a maximum 
at the distance Xfn from the support B, this distance 
being given by Eq. (15). 

Case C. A cantilever beam (Fig. 5) is imiformly 
loaded with a distributed load q per unit length. Here 
the second derivative which enters in Eq. (1) is posi- 
tive, hence if the expression for the moment be also 
considered as positive no mmus sign is needed m the 
equation. At a distance x from the wall, or Xi from the 
free end, the bending moment is 

^ = ?^i-|=i?(?-^)S . . . (17) 

so that Eq. (1) becomes 

g=2i/('-*- • • • • ('« 
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Integrating this expression twice, we get 

dx 2£/L 3 J' • ■ ■ • (19) 

j[Cx + £^' + C,]. . . (20) 



2£/L 



Fio. 5. — A cantilever beam uniformly loaded. 

The constants of integration, C and Ci, are deter- 
mined from the condition that when a; = 0, y = 0, 
and dy/dx = 0. Hence 

C - i P = 0, and i^ i* + Ci = 0; 
or 

The maximum deflection is evidently at the free end 
of the beam, and Eq. (20) gives 

.-C. . . (21) 



18 
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Case D. Let us now consider a case similar to Case 
By except that the left-hand end of the beam (Fig. 6) 
is rigidly fixed in a wall, and the right-hand end rests 
on a support which is lower than the left-hand sup- 
port by a given small amount 5. It is shown in the 
lower part of Fig. 6 that the wall may be replaced 




Fig. 6. — A beam fixed at one end and freely supported at 

the other end. 

by a knife-edge support, provided that the beam be 
extended by a length m and loaded with a weight R. 
The values of R and m must be so chosen that the 
elastic line has a horizontal tangent at A. If this 
condition is fulfilled the two beams in Fig. 6 are sub- 
jected to the same stresses and their elastic curves are 
identical. 
Let the imknown reactions of the supports be Q and 
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Qi. Then the conditions for equilibrium of the forces 
appUed to the beam, and for equilibrium of their 
moments, are: 

P + R = Q + Qi] 

R{m + l) +Pc== Qil; .... (22) 
Pci = Rm + Ql. 

The second of these equations expresses the condi- 
tion that the moment of forces R and P with respect 
to point By which moment tends to rotate the beam 
clockwise, is equal to the moment of force Qi which 
tends to rotate the beam counter-clockwise. The 
third equation expresses a similar condition with re- 
gard to point A. Eq. (6) for the right-hand side of the 
beam becomes 

-^/g = Qx, .... (23) 
and being integrated twice in succession gives 

and 

Ely = Cx - I Qx« + hEI. . (25) 



El'f = C-hQx\ .... (24) 
ax 



The second constant of integration is equal to hEI, 
because y must be equal to 6 when x = 0. 
For the left-hand side of the beam we have 

-^/g = Qit-R{m-\-t), . . . (26) 

which, being integrated, gives 

EI^^=R{mi+\t^)-\Q,t\ . . (27) 

and 

Ely = ft (i w<« + i fi) -I Qx^ . . (28) 
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Both constants of integration are equal to zero, be- 
cause for < = 0, 2/ = and dy/dt = 0. 
Applying now Eqs. (11a) we get 

Eqs. (22) and (29) contain in five equations five un- 
known quantities, viz., Q, Qi, R, m, and C, which can 
thus be detennined. Then Eqs. (25) and (28) will 
give the shape of the elastic curve, and finally the 
abscissa corresponding to the maximum deflection can 
be determined from Eq. (24) by equating its right- 
hand member to zero. This gives 



a:« = V2 C/Q. 

It will be noted that in the case just considered the 
strength of the beam (see Chapter I) can be determined 
only in connection with its elastic curve, because 
Eqs. (22) are not sufficient to determine the reactions 
of the supports. Eqs. (29) contain additional data 
for the restraining moment Rm, for the reactions Q 
and Qi, and for the constant C which determines the 
shape of the elastic cm^e. 

The bending moment M varies along the beam as is 
shown by the broken line in the upper part of Fig. 6. 
It has a positive ma.ximum at the point of appUcation 
of the load, and a negative maximmn at the left-hand 
support. We have 

M,^c = Qc, j . . . . (30) 
Mt^o = -Rm.) 

Either of these moments may be larger in absolute value 
than the other, according to the position of the load 
P; the dangerous cross section may be either at P, or 
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at A. The elastic curve has a point of inflection at a 
point for which the bending moment is zero. Accord- 
ing to Eq. (26), the abscissa U of this point is deter- 
mined from the condition 

Qih = R{m-\-U) (31) 

To the right of this point the beam is concave on the 
upper side, and to the left it is concave on the lower 
side. 

Prob. I. — Deduce expressions corresponding to Eq. (5), for the 
cases in which the load is uniformly distributed, (a) over the whole 
span, and (6) over a portion of the span. 

Prob, 2. — For a beam freely supported on both ends draw to 
scale an elastic curve produced by a imiformly distributed load. 
On the same curve sheet draw an elastic curve for the same beam 
when a concentrated load is appHed at the center of the span. 
Select the loads in such a way as to have the same maximmn deflec- 
tion in both cases. Use a much larger scale for ordinates than for 
abscisssB. Explain from a physical point of view, the reason for 
the difference in the shape of the two curves. 

Prob. 3. — Deduce an expression for the elastic curve when two 
concentrated loads are appHed to the beam, as in Fig. 1. 

Prob. 4. — Modify Case B by adding a load uniformly distributed 
over the span. 

Prob. 5. — Modify Case B by adding a load uniformly distributed 
over a portion of the span. Begin with the simplest case in which 
the distributed load extends from one support to the point of ap- 
plication of the concentrated load. 

Prob. 6. — Deduce an expression for i/max for a cantilever beam 
loaded with a concentrated load at the free end. 

Prob. 7. — Extend Case C by adding a concentrated load at a 
distance c from the wall. 

Prob. 8. — Modify Case C by assuming the beam to be fixed in 
the wall not horizontally, but at a small angle a upward. 

Prob. 9. — A concentrated load travels along a cantilever beam. 
Deduce an expression for the deflection at the end of a beam as a 
function of the position of the load. . 



22 



ENGINEERING MATHEMATICS. 



[Chap. II. 



Prob. 10. — Make all the possible simplifying assumptions in 
Case D, and write, in their final fonn, formulae for the elastic curve 
and for the reactions of the supports. Hint, — Take 6 = 0, and 
take the load at the center of the span. 

Prob. II. — In Case D, assuming 5 to be equal to zero, find a 
position of the load such that the reactions of the two supports will 
be equal. 

Prob. 12. — A horizontal beam is fixed at both ends and is loaded 
with a concentrated weight at the center of the span. Find the 
equation of the elastic curve. 

Prob. 13. — Solve Prob. 12 when the load is not at the center 
of the span. 

Prob. 14. — Solve Prob. 12 when the load is uniformly distrib- 
uted, (a) over the whole span, (6) between a support and a given 
point in the span. 

Prob. 15. — In Probs. 11 to 14 find the position of the point of 
inflection. 

Prob. 16. — In Fig. 6 let 6 be zero when the beam is not 
loaded, and let the support B be elastic, so that when the beam is 
loaded 5 is proportional to Q. Deduce the equation of the elastic 
curve, and write the equations from which the constants are to be 
determined. 




Fig. 7. — A beam with projecting ends, uniformly loaded. 

Prob. 17. — A beam is resting on two supports (Fig. 7), and 
extends over each support by an amount m. It is uniformly 
loaded over its whole length: find the equation of the elastic curve. 

Prob. 18. — For the beam in Prob. 17 find the points of inflec- 
tion and the points where the moment is a maximmn. 

Prob. 19. — Solve Prob. 17 when the load, instead of being dis- 
tributed, is concentrated at C and 2), and at some intermediate 
point between A and B. Take all the three loads different. 



CHAPTER III. 
FATIGUE OF MATERIALS. 

Let a steel bar or rod be placed in a testing machine 
and subjected to a slowly and gradually increasing 
tension imtil it breaks. Let 30 tons be required to 
ruptin-e the bar; this figure is called its ultimate static 
strength. The same bar may be broken by a tension 
of only about 15 tons, provided that the force is ap- 
pUed a great many times and that it varies between 
zero and 15 tons. With such a fluctuating load, pecul- 
iar molecular changes take place in the metal which 
reduce its strength. This phenomenon is known as 
fatigue in the metal. The bar may be broken by a 
load of only 10 tons, if this load varies a great many 
times in succession between 10 tons tension and 10 
tons compression, in other words, between +10 tons 
and —10 tons. 

The effect of fatigue on the ultimate strength of 
steel is represented graphically m Fig. 8. The ordi- 
nates of the curve ABC represent the breaking force 
P as a function of the range of fluctuations of the 
load; the abscissse give the lower limit p of the fluc- 
tuations of the appUed force. Thus, at A, which 
corresponds to, say, 30* tons, the load is stationary, 
p = P, or Oa = aA, and point A is foimd as inter- 
section of the curve with the straight line OA in- 
clined at an angle of 45 degrees to the axis of abscissae. 
At an intermediate point, such as Z), the bar breaks 
imder a load rfZ), when it fluctuates between Od and 
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dD. At point B the load fluctuates between zero and 
OB; this corresponds to 15 tons in our example. At 
G the bar is subjected to alternate tension gG and 
compression Og, the total range of variation of stress 




Compifiasioa Q 



Or Tensioa 



Fig. 8. — Curve of ultimate strength and the effect of fatigue. 

being Gg + Og. By drawing line OC at an angle of 
45 degrees to the axis of abscisssB the point C of the 
curve is foimd, for which the tension Cc and the com- 
pression Oc are equal; in our example cC = Oc = 10 
tons. Parts of the curve outside the range CA do not 
enter into the discussion in this chapter. 

To know the shape of the curve ABC for a given 
material is of the utmost importance in the design 
of machine parts or of a bridge subject to variable 
stresses and shocks. Experimental data are compara- 
tively few and diflScult to get; moreover they cover 
only a few points on the curve; hence the numerous 
attempts of various investigators to estabUsh an em- 
pirical equation of the curve ABC. The purpose of 
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this chapter is to show the reasoning which leads to 
these formulae; this topic furnishes an excellent exam- 
ple of approximate empirical relations established for 
use in engineering practice, where exact data or theory 
are lacking. 

Launhardf s Formula. 

Lamihardt's formula covers only part AB of the 
curve, where the maximiun stress P and the minimunx 
stress p are both tension or both compression. Let 
the ordinate aA be denoted by u, and the ordinate OB 
by V. Launhardt makes a reasonable assmnption that 
the increase of P over v is directly proportional to the 
ratio P : p, or 

P-v = k^, (1) 

where A; is a constant. For the point A, P =^ p = % 
so that the above equation gives for the value of k 

u — V = k. 
Substituting this value of A; in Eq. (1), we get 

P==v + {u-v)^ (2) 

This is Launhardt's formula for stresses of the same 
kind. For practical purposes it is convenient to leave 
the formula in this form, and not solve it for either 
p or P, because the ratio p : P is known before the 
stresses themselves are determined. As the value of 
P represents the breaking strength, to obtain the per- 
missible working stress, P must be divided by a factor 
of safety, which varies from 3 to 10, according to the 
importance of the part, the material used, and the 
certainty with which the stresses can be determined. 
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Weyrauch's Formula. 

By a reasoning similar to that, given above, Wey- 
rauch deduced a formula applicable to the part BC of 
the curve, in other words, when the stress varies be- 
tween a tension P and a compression p. Let the 
ordinate cC be denoted by w. Weyrauch assimies 
that the decrease in the value of P below that of t; is 
proportional to the value of the ratio p : P, or, analo- 
gously to Eq. (2), 

P = e;-(t;-t£;)|, .... (3) 

wherein p is imderstood to represent the absolute 
value of the stress, so that the ratio p : P is always 
positive and less than unity.* For the point J5, p = 
and P = v] for the point C, p =^ P and P = w. For 
other points, P has values intermediate between v 
and w. 

Johnson's Formula. 

An objection to the foregoing two formulae is that 
each of them covers only a part of the useful range of 
the curve ABC, so that there is a break in the slope of 
the ciu-ve at B; in fact, the cm^e ABC is supposed to 
consist of two separate curves, merely having a com- 
mon ordinate v at B. It is more reasonable to assmne 
that both branches are parts of the same curve; con- 
sequently, several formulae have been proposed to cover 
the whole range ABC. The late Prof. J. B. Johnson 
assumed that, in view of inadequate knowledge con- 
cerning the values of v and w, it is suflBicient at present 

* K it is preferred to consider P negative when it represents com- 
pression, Wejrrauch's formula must be written in the form P == v + 
(v — w) p/P. 
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to replace the curve ABC by a straight Ime. This 
line must pass through A, since the value of t^ is known 
with considerable certainty. The slope of the line 
must be such as to satisfy as far as possible the known 
values of v and w. 

The general equation of a straight line between P 
and pis 

P = ap + 6 (4) 

As in the two preceding formulse, it is desirable to ex- 
press P as a function of the ratio p : P. Dividing both 
sides of Eq. (4) by P, we have 



or, solving for P, 






The coefficient 6 can be eliminated because of the con- 
dition that P = u when p : P = 1. Thus, applying 
Eq. (5) to point A, we get 



u = 



1 -a 
hence 

P - <'-")'' (6) 

This is Johnson's formula in its most general form. 
The value of a must be determined experimentally. For 
steel, so far as is known, v is equal to about | u ; the 
substitution of this value in Eq. (6), when p = 0, gives 

I tt = (1 — a) w, 
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or a = i. Hence, for steel, 

p = ^ (7) 

This is the form in which Johnson's formula is usually 
given. This formula gives a value for w oi ^u, be- 
cause P = ^u when p : P = — 1. 

Unwinds Formula. 

Unwin assiunes that the curve ABC forms part of a 
parabola of the second degree between p and P. The 
general equation of such a parabola may be written in 
the form 

{mP + npy +u{rP + sp) +u^f = 0, . . (8) 

where m, n, r, s, and / are coefficients to be determined 
from experimental results. The second and the third 
terms are multiplied by known quantities u and u^, 
respectively, in order to, preserve the homogeneity 
of the equation and to make the coefficients m, n, r, s, 
and/ abstract numbers (numerics). 

Instead of assuming that P depends on the ratio 
p : P, as in the preceding formulae, Unwin assmnes 
that P depends on the range P — p of the fluctuation 
in stress. Let P — p be denoted by A. Then 

mP +np = mP + nP — nP + np = (m + n) P — nA ; 
rP + sp = rP + sP - sP + sp = {r + s)P-' sA. 

Eq. (8) takes the form 
(m+nyP^+n^A^-2n{m+n)PA+u{r+s)P-usA+u^f=0. 

This is an ordinary quadratic equation with respect to 
P. Being solved it gives 

P = QA+Ru + S Vu^ - VuAy . . (9) 
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where Q, R, S, and V are certain combinations of the 
coefficients m, n, r, s, and /. 

If any four points on the curve ABC are known from 
tests, the coefficients Q, R, S, and V may be determined 
by substituting the known values of P and A in Eq. (9) 
and solving the foiu* equations for these coefficients. 
According to the information available at present the 
results of many experiments check satisfactorily with 

formula 

P = i A + Vu^ - VuAy . . . (10) 

in other words, in formula (8), Q = h R = 0, S = 1. 
The value of F is about 1.5 for ductile iron and steel. 
For hard materials, Unwin advises the use of 2 for V. 
With more knowledge of the properties of materials, 
the values of Q, JB, S, and V in Eq. (9) will be determined 
with greater precision. 

Other Possible Formuls. 

For several reasons, it does not seem advisable to 
use the complicated formula (9) to express the curve 
ABC. The formula is imsynametrical, and the ulti- 
mate strength P is expressed through the range of 
fluctuations P — p = A, while it is much more conven- 
ient for practical purposes to express it through the 
ratio of fluctuations p : P. A wide variety of possible 
shapes of the curve ABC may be covered by slightly 
modifying the simple formulae (1), (3), or (5). Thus, 
if experience should show that Eq. (1) does not repre- 
sent the behavior of a metal with sufficient accuracy, 
a quadratic term may be added on the right-hand side, 
so that the formula becomes 

P-v = h(^) + h(^J. . . . (11) 
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With a proper selection of the coefficients A^i and ^2, 
ahnost any degree of curvature may be obtained. 
Weyrauch's formula may be modified in the same way. 
It is also possible to combine Launhardt's and Wey- 
rauch's formulae by using a series of the form 

P = A+B(^) + c(|)'+etc. . . (12) 

A few such terms will represent the results of experi- 
ments with sufficient accuracy. 

Johnson's formula (5) may be made more flexible 
and adaptable in several ways, for instance by adding 
an extra term in the numerator, or in the denomina- 
tor, or in both. These, being merely correction terms, 
need not be calculated with extreme precision. We 
thus obtain the following empirical formulae: 

P ^ (13) 

1-«(|) 

Various other modifications may be introduced by 
using higher powers of p : P. The coefficients which 
occur in these formulae must be determined by using 
a. few known experimental values of P. Thus, for 
instance, let it be required to determine the coeffi- 
<5ients a, b and c in formula (13), when the only ex- 
perimental data available are the values of u, v and w. 
It is convenient in this case to express v and w as 
Jmown fractions of tt, thus 

V = mu; w = nu (15) 
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Substituting these values in succession in Eq. (13) we 
obtain three equations 



u 


"1-a' 


mu 


= 6; 




h-c 



nu = 



1 +a 



(16) 



These equations can be easily solved for a, 6 and c, 
and the results used in Eq. (13) for obtaining inter- 
mediate values of P. Eqs. (11), (12) and (14) may be 
treated in a similar manner. If the nmnber of known 
points on the curve is larger than the number of coeflBi- 
cients in the empirical formula, the most probable values 
of the coefficients may be determined by the method 
of least squares or by trial. See Part III (Thermody- 
namics), Chapter III. 

Prob. I. — Plot values of P according to Launhardt's and Wey- 
rauch's formulae for v = iu and w = \u. Also draw the straight 
Kne corresponding to Johnson's formula, and determine graphically 
maximum per cent discrepancy between the curves and the straight 
Une. 

Prob. 3. — Check the above maximmn per cent discrepancy 
analytically, by equating to zero the first derivative of the analytical 
expression for per cent discrepancy. 

Prob. 3. — There is a discontinuity in the slope where Laun- 
hardt's curve passes over to Weyrauch's curve. Determine ana- 
lytically the angle between the two curves. 

Prob. 4. — Show that if w + lo = 2 v, Launhardt's and Wey- 
rauch's formula? represent parts of the same ciure. 

Prob. 5. — A certain kind of steel breaks at a steady tension of 
40 kg. per sq. mm. In order to break a similar sample by a tension 
of 20 kg. per sq. mm., tension must alternate with a compression of 
10 kg. per sq. mm. Write Johnson's formula for this material. 

Prob. 6. — Plot curves ABC according to Unwin's formula (10) 
for ductile steel and for hard steel. Draw straight lines that agree 
as nearly as possible with these curves and calculate the correspond- 
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ing valuer of a in formula (6). Assume u to have the same value 
in all four cases, the four curves passing through the same point A. 

Prob. 7. — For a certain material v =« 0.62 u and w = 0.45 u. 
Can the properties of the material be satisfactorily represented by 
formula (10), and if not, what will be the numerical coefficients in 
the more general formula (9)? 

Prob. 8. — Using relation (16) determine the values of the coeffi- 
cients in formula (12), in terms of m, n and u, 

Prob. 9. — Using relation (15) determine the values of the coeffi- 
cients in formula (13), in terms of m, n and u, 

Prob. 10. — Using relation (15) determine the values of the coeffi- 
cients in formula (14), in terms of m, n and u, 

Prob. II. — From the available literature on the subject, plot 
experimental fatigue curves, such as ABC^ for various materials, 
and investigate which of the formula quoted in this chapter can be 
made to express the results of the tests most accurately. 



CHAPTER IV. 

EULER^S THEORY OF FLEXURE OF LONG 

RODS. 

A LONG elastic rod or flat spring (Fig. 9) is pivoted at 
both ends and is loaded with a force P in such a way 
that it is bent out of its original straight shape. A 
guide g constrains the upper end of the rod to remain 
on the original vertical line. The problem is to find 
a relation between the force P and the deflection a 
which it produces. According to the general theory 
of flexure, which the student can find in any standard 
work on mechanics of materials, we have, for any cross 
section of the rod such as A, 

^ = M. ...... (1) 

p 

In this expression M is the bendmg moment of the 
external forces with respect to the cross section under 
consideration, in kg.-cm.; E is the modulus of elastic- 
ity of the material, in kg. per sq. cm. ; / is the moment 
of inertia of the cross section with respect to the axis 
perpendicular to the plane of Fig. 9 and passmg 
through the center of gravity of the cross section, in 
cm.^; p is the radius of cm-vature of the axis of the rod 
at Aj in cm. 

When the curvature is small p is large and we have 
approximately (see Prob. 10 below) 

1 ^ _d^ 

p dx^' 

33 



34 



ENGINEERING MATHEMATICS. 



[Chap. IV. 



1. 



Iti 
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Fio. 9. — A long rod subjected to endwise compression. 



SO that Eq. (1) becomes 



-TE^ = 



dx* 



M. 



In the case under consideration the bending moment 
M = Py, so that 



(2) 
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or, denoting, for the sake of abbreviation, 



V/,^ = "' (3) 



we have 



dx^=-"^V (4) 



This is the differential equation of flexure, and it shows 
that y is such a function of x that its second derivative 
with respect to x is proportional to the function itself, 
the coefficient of proportionaUty — n^ being negative. 
A sine (or cosine) function of the general form 

2/ = a sin (6a; + c) (5) 

is the only function which satisfies this condition.* 
In this formula, a, 6, and c are certain constants to be 
determined from the given conditions of the problem. 
Eq. (5) means that the axis of the rod is bent into a 
sine ciu^e. To find its exact shape we differentiate 
Eq. (5) twice: 

-i^ = ab cos (bx +c) (6) 

dx 



dx^ 



a¥sm{bx +c). . . . (7) 



Substituting the values from Eqs. (5) and (7) into Eq. 
(4) and cancelling the sine term, we obtain 

— ab^ = —n^a, 
or 




^=^=\lE •• • (8) 

* If the right-hand side of Eq. (4) had a plus instead of a minus 
sign, the simplest solution would be of the form y =/e-**. For a 
complete integration of Eq. (4) see the Note at the end of the chapter. 
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The constant c is determined from the terminal condi- 
tions of the problem; namely, when x = or x = I, 
y must be equal to zero. Hence 

sin c = and sin (bl + c) =0, 

or, in the simplest case, c = and hi = x. This shows 
that the rod forms a complete half-wave of a sinusoidal 
curve. Using the value of 6 from Eq. (8), we get 

V§ = ^> (8a) 

or 



t'EI 



P--^ (9) 

The amplitude a of the sine-wave in Eq. (5) remains 
indefinite, because Eq. (2) is satisfied with the solution 
for any value of a whatsoever, and there are no other 
physical or mathematical conditions from which a 
could be determined. This result may be taken to 
indicate that, as soon as the force P reaches the 
critical value, Eq. (9), the rod may bend to any value 
of deflection a, once it has started to bend due to 
the slightest lack of symmetry. One must remember, 
however, that a must be small compared to Z, because 
of the approximate value for p used in Eq. (2), which 
is true for small curvatures only.* 

For practical piuposes Eq. (9) determines the break- 
ing load P, because no part of a structure can be con- 
sidered safe, in which such an indefinite buckling may 
take place. It is convenient to know the crushing force 

* For a 'differential equation of flexure of columns more accurate 
than Eq. (4) see R. F. Deimel, The Theory of Columns^ Bulletin of the 
Society for the Promotion of Engineering Education, Vol. 6, p. 258, 
December 1915. It is shown there that for each load there is a definite 
deflection. This is also confirmed by experimental investigations. 
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per unit area of cross section of the rod. Designating 
this area by A, and dividing both sides of Eq. (9) by A, 
we have 

A P \AJ 

But the moment of inertia /, being of the dimension 
(length)^, may be represented as 

/ = Ar\ 

where the fictitious length r is called the radius of gyra- 
tion of the cross section ; thus, we finally get the well- 
known theoretical formula for the flexure of long rods 

l = J^ (10) 

In the preceding theory the rod is supposed to be 
hinged at both ends. Three other cases of flexure are 
shown in Figs. 10, 11, and 12. The first figure shows 
a rod fixed at the lower end B, and loaded at the free 
end A. This case could be considered in a way similar 
to that outlined above, only the termmal conditions for 
the determination of the constants of integration are 
different. However, the result may be foreseen directly 
by imagining the rod extended dowward, as shown by 
the dotted line. The case is reduced to that in Fig. 9 
and formula (9) is applicable, if I is imderstood to be 
the total length of the imaginary rod between A and C; 
or I = 2 L, where L is the actual length of the rod. 

In a similar manner, the rod in Fig. 11, fixed at both 
ends, B and H, may be extended both ways, by one- 
quarter of the sine-wave, so that we have three rods, 
CDf DE, and EF, one on top of another, bent accord- 
ing to Fig. 9. Eq. (9) is again applicable if the theo- 
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Fig. 10. — A long vertical rod fastened at the lower end and loaded 

at the upper end. 
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Fig. 11. — A long rod fastened at both ends and subjected to a 

vertical load. 
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retical length I be taken equal to \ L, where L is the 
actual length of the rod. 

The f oregomg theory is applicable only to very long 
thin rods, for which the direct compressive action of 
the load P is negligible as compared with the bend- 
ing moment produced by a deflection. For ordinary 
shorter columns and for struts used in bridge construc- 
tion, various semi-empirical formulsB are used. The 
elastic curve of a short colunm is not a sine wave, so 
that the relations between L and I deduced above do 
not exactly hold. The strength of such colxunns is 
discussed in the next chapter. 

Note, — To find the general solution of Ekj. (4), multiply both 
sides of it by dy/dx. This gives 

scen—iw <"> 

or 

(|[y= -nV + nV (12) 

The tenn nV represents the constant of integration, a being so far 
an arbitrary quantity. The constant could be simply written C, 
but in physical problems it is of advantage to give all constants a 
certain physical meaning. Thus, y being a length, a has also the 
dimension of a length. When dy/dx = 0, y = a, so that a is the 
TTiR Yi'TniiTn deflection of the rod. We have further 

^ = ±nV^^^:7 (13) 

ax 

The double sign means iti this case that the rod may buckle either 
to the positive or to the negative side of the vertical axis. Limit- 
ing our attention to the case shown in Fig. 9 we need retain only 
the plus sign. Hence 

or, after integration, 

x-irsin-i^-c], (14) 

nL a J 
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where c is a new constant of integration. Solving this equation for 
y, expression (5) is obtained; it will be remembered that b is equal 
to n, according to Eq. (8). 



Prob. I. — Assuming the crushing 
load in the case indicated on Fig. 9 
to be 100 per cent, what are the values 
of crushing load for identical rods of 
the same actual length L, but fastened 
as shown in Figs. 10 and 11? 

Prob. 2. — A spring of length L is 
bent as shown in Fig. 12, forming N 
sine waves (six half-waves in the par- 
ticular case shown in the figure). How 
many times larger force P is required 
to hold the spring in this position, as 
compared with the force necessary to 
hold the same spring in the shape 
shown in Fig. 9? 

Prob. 3. — For round rods formula 
(9) may be represented in the form 



P = Const 



■(!)■ 



<^ 



P 



H 



where d is the diameter of the cross 
section. Find the numerical value of 
the constant for steel rods, d and I 
being in cm. and P in kg. The mod- 
ulus of elasticity E for steel is about 
2 million kg. per sq. cm. 

Prob. 4. — Compare the crushing 
load on a round rod to that on a Fig. 12.— A long spring held 



m 



7/ 



by an end force so as to 
form six half-waves. 



rod having a square cross section of 

the same area. If the prescribed 

crushing load be the same in both cases, which rod is lighter and 

by how much? 

Prob. 5. — A flat rod of cross section 0X6 (Fig. 13) is to be 
strengthened by two ribs of thickness c, so as to be as rigid against 
buckling in the direction XX as it is in the direction YY. Deter- 
mine the height z of the ribs. 
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Prob. 6. — Deduce the expressiona for the radius of gyration of 
some simple croaa sections used in practice, such as are shown in 
F^. 14. 



Pio. 13. — A rectangular eroes section atreogthened by ribs. 
Prob. 7. — Sketch the general shape of the curve which gives 
the unit stresB P/A as a function of l/r, Wliy does the curve seem 




Fig. 14. — Cross 1- 



!■■ -I iins used it 



a practice 

to indicate that for small values of I the rod can apparently support 
an infinite load, which is of course not the case? 

Prob. 8. — Deduce an expression sundar to Eq (9) for the 
case shown in Fig. 10; use a differential equation similar to 
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Eq. (2), and apply the proper terminal conditions for a; = and 
X = L. The result must check with that anticipated in the text 
above by reducing this case to that shown in Fig. 9. Hint. — The 
bending moment is P (a + t/), and one of the terminal conditions 
is that dy/dx = when x = 0. 

Prob. 9. — Deduce an expression similar to Eq. (9) for the case 
shown in Fig. 11; use a differential equation similar to Eq. (2), 
and apply the proper terminal conditions at the ends. The result 
must check with that obtained directly from Eq. (9) by putting 
the effective length I equal to one-half the actual length L of the rod. 
Hint. — For a point on the rod, such as A, the bending moment 
is equal to Py — Jlf 0, where Mq is the so-called moment of fixture. 
This is the moment exerted at the upper end of the rod to keep 
it vertical. The moment of the force P tends to turn the upper 
part of the rod counter-clockwise, while the moment Mq tends to 
turn it clockwise; for this reason Mo enters the equation with the 
minus sign, if Py enters with the plus sign. The terminal conditions 
are: at a; = 0, y = and dy/dx = 0; at x = Lj dy/dx = 0. 
There are, therefore, three equations from which Mq and the two 
constants of integration may be evaluated. 

Prob. 10. — Deduce an expression similar to Eq. (9) for the 
case when the upper end of the rod is hinged as in Fig. 9, and the 
lower end is fastened as in Fig. 11 ; use a differential equation simi- 
lar to Eq. (2), and apply the proper terminal conditions at the ends. 
Consult the two preceding problems. 

Prob. II. — Prove that when the curvature of the rod is small 
(Fig. 9), it is approximately equal to d^/dx^. First proof. — By 
definition, the curvature 1/p = da/dsy where a is the angle which 
the tangent to the curve at a point makes with the X axis, and ds 
is an element of the curve. In general, tan a = dy/dx, but in the 
case imder consideration the angle is small throughout so that ap- 
proximately a = dy/dx. Besides, ds is practically equal to dx, 
80 that 1/p = d (dy/dx) /dx = d^/dx\ Second proof. — Take the 
general expression for the radius of curvature given in the differ- 
ential calculus, and let dy/dx = 0. 

Prob. 12. — In addition to the end load P (Fig. 9) the rod is 
subjected to a uniformly distributed horizontal load of p kg. per cm. 
length. Write down the differential equation for the bending 
moment and integrate it. Solviion. — The total horizontal load 
is plj and its reaction at each end of the rod is i pi. The moment of 
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the upper reaction with respect to point A is equal to i pZ (i — x). 
The load itself over the upp)er portion of the rod is p (Z — x)j and 
its center is distant by § (Z — x) from A. Thus, the moment of 
the load is — i p (Z — xf. Therefore, the right-hand member of 
Eq. (2) must be supplemented by the term 

ipZ (Z - x) - Jp (Z - x)2 = ipx (Z - x), 
Eq. (2) is then modified to one of the general form 

^, + nV=-Ax(Z-x), (15) 

where n and h are known coefficients. According to the general 
theory of linear differential equations with constant coefficients we 
find that 

2/ = (Asinnx + JJcosnx) i(Zx — x*-h~), . (16) 

where the first expression in the parentheses on the right-hand side 
is the " complementary function," and the second expression is a 
'^ particular solution." A and B are the constants of integration, 
and are determined from the terminal conditions, as before. 

Prob. 13. — Prove by a direct differentiation that solution (16) 
satisfies Eq. (15). 

Prob. 14. — With reference to Eq. (8a), show how the value for 
P is obtained corresponding to Fig. 12, by using mw in place of t, 
where m is an integer. 



CHAPTER V. 
FLEXURE OF COLUMNS. 

Euler's theory of flexure given in the preceding 
chapter is only appUcable to rods, colunms, springs, 
etc., whose length is several hundred times greater 
than the transverse dimensions (radius of gyration). 
The reason is that in Euler's theory the direct com- 
pression is left enth-ely out of consideration. This 
compression becomes more and more important as 
the rod becomes stiffer and shorter, that is, as its 
length becomes smaller in proportion to its transverse 
dimensions. 

For short blocks the carrying capacity, indicated by 
formula (9) in the preceding chapter, tends to become 
infinite, while in reality it is evidently limited by the 
ultimate compressive strength P/A of the material. 
Various empirical and semi-empuical formulae are used 
by engineers for calculating the strength of columns of 
medium length, which are too short to be computed 
according to Euler's formula, and at the same time are 
too long to be considered simply as blocks subjected 
to a uniform compression. Some of these formulae 
and their deductions are given below, as examples of 
reasoning by which practical problems are solved 
approximately, where the true physical relations are 
too comphcated. The student will also get some 
experience in fitting empirical formulae to experimental 
curves. 
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(1) Johnson^ 8 Straight-line Formula. — Euler's theory- 
leads to the following relation between the crushing 
load per unit area P/A, and the ratio l/r: 

A {l/ry 



(1) 



See formula (10) in the preceding chapter. It will 
be seen that the crushing stress depends not upon the 




Fig. 15. — Euler^s theoretical curve and Johnson's straight-line 

correction. 

absolute length of the rod or column, but upon the 
ratio of the length to the radius of gyration of the cross 
section. Plotting the values of P/A against '/r as ab- 
scissa a curve is obtained of the general shape shown 
in Fig. 15. This curve gives absurd values of P/A 
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for small values of l/r, as is explained above. Thomas 
H. Johnson proposed therefore to replace the upper 
part ah of the curve by a straight Une ac tangent to it. 
This straight line is determined by the condition that 
for very short rods, that is, when I = 0, the imit stress 
must be equal to the ultimate crushing strength p of 
the material. 

The equation of this straight Une is found analytically 
as follows : Let, for the sake of abbreviation, 

I P 

r A 



2/ = — » (2) 



We have then for Euler's curve 

X' 

and for the straight line 

y = p -kxy (3) 

where A; is a coefficient to be determined from the con- 
dition that the straight line is tangent to the curve (2). 
Let the coordinates of the point of tangency be Xi 
and j/i. Then, at this point the values of y from 
Eqs. (2) and (3) must be equal to each other as must 
also the values of dy/dx. Hence 

IT Hj 

— ■ = p - kxi, (4) 

X\ 

and 

--^=-^ (^> 

From these two equations k and Xi can be determined. 
Substituting the value of k from Eq. (5) into Eq. (4) 
gives 



I 
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or 

-a-v'f («) 

Using this value of Xi in Eq. (5) we get 



^=l^v/ 



and Johnson's formula becomes 

This formula should be used only between the values 
l/r = 0, and l/r determined by the relation (6). For 
longer columns Euler's formula (1) is vaUd. 

(2) HodgkinsorCs Formula. — Hodgkinson was one 
of the first experimenters on the strength of long rods, 
and he foimd that the crushing load does not vary 
exactly as required by Euler's formula (1). He there- 
fore modified it into an empirical formula 

P r r^ 



= C.W, (9) 



where the coeflicients C, m, and n are to be so selected 
for each material, or for each kind of colunms, as to 
satisfy the results of available experiments. Formula 
(9) is not now used, since the experimental data on 
colunms can be better expressed by more rational 
formulae given below. But this old formula deserves 
some attention, because in quite a niunber of engineer- 
ing problems a satisfactory agreement with experi- 
mental results is obtained by the use of fractional 
exponents, in place of integers required by a simplified 
theory. 

(3) Rankine^s Formula. — Side by side with empir- 
ical formulae, several prominent engineers, among 
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others Navier, Rankine, and Ritter, worked on a rational 
theory of flexure of short columns. The theory now 
accepted is as follows: A coliunn of moderate length 
fails mider the combined stress caused by direct com- 
pression and by flexure. The uniform stress due to the 
direct compression is P/A, while the maximum com- 
pressive unit stress pi caused by flexure on the concave 
side (Fig. 2) is 

Jh^-Y (10) 

In this formula M and / have the same meaning as in 
the preceding chapter, and z is the distance from the 
center of gravity of a cross section to the extreme com- 
pressed fiber.* The maximum bending moment M 
is proportional to the maximum deflection a of the 
colmnn (Fig. 9), or 

M = Pa. 

Thus the total compressive stress in the extreme fibers 

on the concave side is 

P , Paz 

or, substituting for I its value Ar^, as in the preceding 
chapter, we finally get 

Consider two similar colmnns, similarly loaded, so that 
their elastic curves are geometrically similar. Then, 
both the deflection and the dimension z of each column 
are proportional to its length, so that 

a = qil, 
z = q2l, 

* For a proof of formula (10) see any standard work on mechanics 
of materials. 
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where gi and 52 are empirical coeflScients. Substitut- 
ing these values in Eq. (11) and denoting the product 
gi?2 by g, we get 

P=j[l+q(l/ryi 



or 

p ^ 



A 1+q (l/r) 



(12) 



This is Rankine's formula, in which g is an empirical 
coefficient and p is the actual maximum compressive 
stress on the concave side of the colmnn. For short 
prisms the second term in the denominator is small as 
compared with imity, so that the formula gradually 
changes into the ordinary formula, P/A = p, for uni- 
form compression. On the other hand, for very long 
columns, imity in the denominator may be neglected 
as compared with the other term, and we get 

A q\l/ry ^-""^^ 

(4) Ritter^s Formula. — In Eq. (12) the coefficient q 
may be so determined that for very long columns this 
formula would give the same value of buckling load as 
Euler's formula (1). In Euler's formula P/A is the 
stress in the material up to or beyond the elastic limit 
Pe, so that equating the right-hand members of Eqs. (1) 
and (13) imder this supposition, we have 

J[!^=l_2^ (14) 



whence 



^-^ (IS) 
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Formula (12) then becomes 

f ^ (16) 

1 + ^ (l/ry 

In this form it is known as Ritter's formula. 

It must be clearly understood that Rankine's for- 
mula (12) may be used for determining either the safe 
load or the crushing load of a coliunn, as long as q 
remains an empirical coefficient, determined directly 
from experiments. However, when q is ascribed a 
theoretical value, as in formula (16), the actual stress 
p cannot exceed the elastic limit of the material, 
as will be seen from the method of deducing the for- 
mula. Therefore, formula (16) cannot be made to check 
directly with results of experiments on failure of col- 
mnns. However, the same objection may be raised 
against many rational formulae used in engineering 
work without in the least impairing their usefulness for 
design. 

(5) VierendeeVs Formula. — Prof. A. Vierendeel of 
the University of Louvain (Belgimn) has shown that 
a formula similar to Rankine's may be deduced di- 
rectly from Euler's formula, by considering the mod- 
ulus of elasticity -EJ to be a variable quantity. Namely, 
experiments show that beyond the elastic limit the co- 
efficient of proportionality E between the stress and 
the strain decreases with the increase in the actual 
stress to which the material is subjected. For steel 
and iron a straight-line law may be assimied 

E = a-h^, (17) 

where a and 6 are empirical coefficients. Vierendeel 
states that he is justified in using this empirical ex- 
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pression for E in formula (1) since Euler's formula is 
supposed to give a critical load, which strains the 
material almost to its ultimate strength. Thus, he 
obtains 

P TcHa-bP/A) 
A {l/rY 



or, solving for P/A, 

P T^a 



A ir% + {l/rY 



(18) 



This formula is mathematically identical with formula 
(12), though the coefficients a and h have an entirely 
different physical meaning. Formula (18) is intended 
for columns of mediiun length, for which the critical 
imit stress P/A comes out beyond the elastic limit of 
the material. For long columns Euler's formula (1) 
must be used. The value of JS? in Eq. (1) and the values 
of a and h in Eq. (15) must be so selected that there is 
no discontinuity in the values of the critical load, and 
both formulae must give the same result when P/A 
is equal to the elastic limit of the material. 

Vierendeel's formula cannot, strictly speaking, be 
considered a rational formula, but it checks quite 
closely with the results of munerous experiments and 
for this reason ought to be useful in practice. An ad- 
vantage of this formula over Rankine's is that it con- 
tains two empirical coefficients, instead of one, and 
therefore can be made to represent the results of experi- 
ments much more acciu'ately. 

Prob. I. — The curve in Fig. 16 gives the results of certain ex- 
periments on the buckling of wrought-iron columns with pivoted 
ends. Show how closely Euler's curve covers the points corre- 
sponding to high values of l/r^ and determine Johnson's straight 
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line (both analytically and graphically) for low values of l/r. Take 
E = 1.75 million kg. per sq, cm., and p =■ 3000 kg. per sq. cm. 

Prob, 3. — Prove that the ordinate of the point of tangency 
between Eider's curve and Johnson's straight line always corre- 
sponds to one-third of the ultimate strength p of the material. 



Ratio I-M- 

Fia. IS. — An eiqierimental curve of critical load for certain wrought- 
iron columns with hinged ends. 

Prob. 3. — Suppose that the straight line ca (Fig. 15) does not 
represent the results of certain experiments with sufficient accu- 
racy, so that a more compUcated curve is required to cover the range 
between a and c. Let this curve be required to pass through a given 
point ft {xd, ^o)- Suggest a few simple possible equations for the 
curve cka, for instance y ^ p — kixT ktx^, or y = 7-— — , and 

show how to determine the values of the coefficients in these 
equations in order that the curve be tai^nt to £uler's curve. 
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Prob. 4. — Represent the experimental curve shown in Fig. 16 
by an equation of the form 

P C 

A (l/r)"^ ' 

select the coefficients C and m so as to cover the widest possible 
range of values of (l/r). Try logarithmic cross section paper. 
Consult the article on the most probable curve in Part III 
(Thermodynamics), Chapter III. 

Prob. 5. — For very long colimins the critical load varies in- 
versely as the square of the length of the column. Show by means 
of the experimental curve and of the formulae given above that 
for columns of medium length the variations in strength are less 
than the square of the length, and that for very short prisms the 
critical load is independent of the length of the block. 

Prob. 6. — Determine p and q in Rankine's formula so as to 
cover as well as possible the values of critical stress for the lower 
and medium values of l/r in Fig. 16. 

Prob. 7. — Draw a curve of load per sq. cm. of cross section that 
would strain the columns in Fig. 16 up to the elastic limit, accord- 
ing to Hitter's formula. The elastic limit of the material is about 
one half of the crushing stress given in Prob. 1. Keep in mind 
that Hitter's formula does not give values of crushing loads, but of 
safe loads only. 

Prob. 8. — Express the experimental results shown in Fig. 16 
by means of Vierendeers formula for medium values of Z/r, and by 
means of Euler's formula for large values of l/r. Find the value of 
l/r beyond which Euler's formula must be used. Are the two 
curves tangent to each other, and if not, what is the angle between 
them? 

Prob. 9. — The external diameter of a hollow circular column is 
38 mm., the internal diameter is 25 nun., and the length is 25 m. 
The ends of the column are provided with flanges and are securely 
bolted to the floor and ceiUng beams. What value of l/r should be 
used in calculating the strength of the colunm? Consult the pre- 
ceding chapter in regard to the effective length of a column with 
the ends fastened. 



CHAPTER VI. 

STRENGTH OF THICK CYLINDERS AND 

SPHERES. 

Let Fig. 17 represent a cross section of a steel 
cylinder with thick walls. Such cylinders are used for 
hydrauUc presses, also as containers for compressed 
gases, and as barrels for big guns. Let the cylinder be 
subjected to a given internal pressure po kg. per sq. 
cm. due to a gas or a Uquid. The problem is to find 
the mechanical stresses in the wall, or to determine 
the thickness of the wall such that the maximum stress 
will not exceed a prescribed safe hmit. 

The stresses within the wall are tangential, that is, 
perpendicular to radii, and oppose the blunting of the 
cylinder. * If the cylinder be thought of as consisting 
of concentric shells of infinitesimal thickness, the stress 
varies from shell to shell; the tangential tension is a 
maximum at the inner surface and reaches its minimum 
at the outer surface of the cylinder. The difficulty in 
the problem is to find the exact law of distribution of 
stress within the wall, so as to be able to calculate the 
stresses. 

This problem of finding the distribution of stress in 
a thick cyUnder is of interest to a student of applied 
mathematics, because he can follow in it a gradual de- 
velopment of the theory, and he can compare the im- 
provements made in the fundamental assumptions. 
Five successive assumptions or approximations are 
Inade in the treatment given below, and the results so 
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derived are analyzed. The first approximation is the 
simplest and the crudest, the last is in accordance with 
the modem views on the theory of elasticity, though it 
leads to much more complicated mathematics. 




Fig, 17. — Ctobb section of a thick hollow cylinder or sphere. 

The problem of the strength of a thick hollow sphere 
subjected to internal pressure is similar to that of the 
cylinder, and is briefly treated in parallel with it, in 
order to give the student an opportunity of applying 
the various assumptions made for the cylinder to a 
body of a different geometric shape. 

Let the thickness of the wall (Fig. X7) be S em., the 
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inner radius ro cm., the outer radius n em., and the 
internal pressure po kg. per sq. em. We shall assume 
at first that this internal pressure is so great that the 
external atmospheric pressure pi may be neglected in 
comparison. 

Should the cylinder burst imder the influence of this 
internal pressure, the two parts of it would separate 
along a plane such as AB, so the condition for equi- 
libriiun is that the siun of the upward vertical com- 
ponents of gas pressm-e on the upper half of the 
cylinder is equal to the siun of the vertical stresses 
along AA^ and BB\ 

Consider a slice of the cylinder one cm. long in the 
direction perpendicular to the plane of the paper. 
The surface corresponding to an infinitesimal central 
angle da is roda, and the normal pressure is po^oda. 
But for each element such as ro da there is a symmetrical 
element on the other side of the plane OCy so that the 
horizontal components of the pressure are balanced 
and the vertical components are added. The vertical 
component of the pressure upon either element is 

Po^o da COS a:, 

SO that the total pressure on either quadrant is 

r 

Po^o cos ada = PoVq. 



f 

Jo 



Since ro X 1 is the area of the plane A'O, the foregoing 
result may be interpreted by saying that the total 
vertical pressure is the same as if it were distributed 
uniformly upon the diametral plane of the cylinder. 

Let the variable unit stress along AA' be denoted by 
t, and the distance from the axis by x. Then the in- 
finitesimal force acting upon an element of the plane 

4 



*/ro 
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AA', per unit of length in the direction perpendicular 
to the plane of the paper, is t dx. Equating the total 
internal force resisting a rupture along AA' to the total 
pressiu-e of gas, we get 

tdx^' Po^o (1) 

To 

This is the fundamental equation for the cylinder. 
In order to be able to integrate the left-hand side it is 
necessary to make an assumption regarding the dis- 
tribution of < as a fimction of the distance x. This is 
done in the various approximations that follow. 

Let now Fig. 17 represent the cross section of a 
thick sphere; the area of the diametral plane A'fi' is 
TTTo^y and the upward pressure of the gas is poitr^?. The 
stress corresponding and perpendicular to an infini- 
tesimal strip of plane A A' of width dx is ^ • 2 xx • dx. 
Thus, equating the pressure of the gas to the internal 
stress produced by it we have 



/ 



t '2tx 'dx = PotTq^, 



f 



or, dividing by 2 t, 

txdx = ^ PoVo^ (2) 

This is the fundamental equation for the sphere. A 
reasonable assumption must be made concerning ^ as a 
fimction of x, before this equation can be integrated. 
Some of such assumptions are discussed below. 

First Approximation. 

The Stress is Distributed Uniformly. 

(a) Cylinder. — The simplest assmnption that can 
be made in Eq. (1) is 

t = to =^ const. ; (3) 
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in other words, the stress is distributed uniformly 
through the thickness of the wall. This assumption 
is practically correct for thin walls, for instance, in 
steam boilers, and it may be accurate enough for 
comparatively thick walls, where particular accuracy 
is not required. 

Substituting a constant value of to in Eq. (1), and 
integrating, we get 

to in - To) = pon, 
or 

to8 = poro, (4) 

where d is the thickness of the wall. Any one of the 
four quantities which enter this formula may be de- 
termined if the remaining three are known. 

If the external pressure pi cannot be neglected, 
Eq. (4) becomes 

to8 = poro - Pin (4a) 

(6) Sphere. — Using condition (3) in Eq. (2) we get 

to (ri2 - ro^) = poroS 
or 

W (n + To) = PoTo^ (5) 

For thin walls n is approximately equal to r©, and the 
formula is simplified to 

2to8 = poro (5a) 

Prob. I. Prove Eq. (4a) and explain the reason why 5, accord- 
ing to this formula, assumes a negative value when the external 
and the internal pressures are equal. 

Prob. 2. Instead of bursting along a diametral plane, a cylinder 
may burst along a cross section perpendicular to its longitudinal 
axis, due to pressures on the enclosed ends. Show that imder the 
assumption of a imiform distribution of stress and thin walls, a 
cylinder possesses a double factor of safety against such a rupture 
as compared to a failure along a diametral plane. Hint. — Let the 
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unit stress parallel to the longitudinal axis be q kg. per sq. cm.; the 
resisting stress is approximately equal to 2irrifiq, The pressure of 
the gas on the ends is xTq^o. 

Prob. 3. Apply Eqs. (5) and (5a) to the case in which the exter- 
nal pressure pi cannot be n^lected. 

Second Approxiniation. 

. The Thickness of the Wall is Assumed Unchanged 

Under Stress. 

(a) Cylinder. — A more correct idea of the dis- 
tribution of the stress t throughout the cylinder wall 
is obtained by analyzing critically the foregoing as- 
sumption of t = constant. According to the funda- 
mental law of the theory of elasticity, stresses are 
proportional to deformations which they produce 
(Hooke's law). Let I be the original length of a fiber 
subjected to a unit stress t, and let AZ be the elongation 
of this fiber under the influence of t. Then the relative 
elongation (or the elongation per unit length) is Al/l, 
and the above law is expressed mathematically thus : 

t = Ef, (6) 

where E is the so-called coefficient (or modulus) of elas- 
ticity of the material. This law holds true as long as 
the elastic limit of the material is not exceeded. Apply- 
ing this law with the assumption that t is constant, 
it is shown below that one gets the absurd result that 
imder the influence of an internal pressure various 
concentric annuli of the cylinder tend to spread apart. 
Common sense would indicate that an internal pres- 
sure tends to press the inner annuU or layers against 
the outer ones, because only in this way can the outer 
layers be made to help in the resistance to bursting. 
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Let u be the increase in the radius x of an infini- 
tesimal shell within the wall (Fig. 17), under the in- 
fluence of the internal pressure; in other words, the 
particles of the material which originally are at a dis- 
tance X from the center of the cylinder occupy a circle 
of radius x + u. According to the above mentioned 
law of elasticity, the tangential stress to which these, 
particles are subjected is 

<=g. 2Wx + «)-2.a; ^ ... (7> 

2 tx 
or 

t = E- (7a) 

X 

If t be assumed constant, as in the first approximation, 
then u must be proportional to x; this means that the 
external annuli expand more than the internal ones, 
all the annuh or layers tending to spread apart, which 
is evidently absurd. In reahty u decreases with the 
increasing x, so that the internal annuli press upon the 
external ones. Not knowing the actual law relating 
u and X, the next best assumption is that u is constant, 

u =c; (8) 

in other words, there is no spreading or compression be- 
tween the concentric layers and the total thickness of the 
wall remains unchanged by the pressure. Eq. (7a) then 
indicates that t varies inversely with x, so that the 
inner parts of the wall are subjected to greater stresses 
than the outer parts. The results thus obtained, while 
not altogether correct, must be more accurate than the 
relationship stated in Eq. (4). 

Substituting the value of t from Eq. (7a) into Eq. (1) 
and integratmg gives 

EuLn- = poro (9) 

^0 
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Eliminating the product uE from Eqs. (7a) and (9) 
we obtain the following formula for t: 

The maximum tensile stress is at the inner surface of 
the cylinder where x = r©, and is equal to 

^°*" Ln(n/ro) ^''^ 

(&) Sphere. — It is assumed in Eq. (7) that all the 
stresses are in planes perpendicular to the longitudinal 
axis of the cylinder, so that the value of the stress 
depends upon the linear expansion of the wall. In the 
case of a sphere the stresses at a point have the same 
value in all directions in a plane perpendicular to the 
radius. Each infinitesimal concentric spherical shell 
remains a sphere after expansion, so that stresses may 
be said to depend upon the surface eTpansion. By 
analogy with Eq. (7), we may write 

where 4 tx^ is the area of the infinitesimal shell within 
the wall. This equation is simplified to 

2ux +u^ 



t = E 



X 



2 



We consider here deformations u which are extremely 
small as compared with the radius x of the shell. 
Therefore u^ in the numerator of the preceding formula 
may be neglected without any appreciable error, and 
we get 

i = E^ (12a) 

X 
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Eq. (12a) may also be written directly from Eq. (7), 
by considering a smaU square piece in the infinites- 
imal shell imder consideration. This element is sub- 
jected to stresses t on all four sides adjacent to the rest 
of the shell, and each pair of equal and opposite forces 
t not only stretches the element, but also coimter- 
balances the contracting action due to the other pair 
of forces t. Assume that it is necessary to apply a 
tension stress equal to ^ < to counterbalance the sec- 
ondary compression eflfect due to a tension stress t in 
a perpendicular direction. The conclusion is that all 
the stresses t must be doubled for a given u in a, 
spherical shell, as compared with a cylindrical shell. 
It will be seen accordingly that t has a double value 
in Eq. (12a) as compared with Eq. (7).* 

Proceeding as in the case of the cylinder we obtain 

Prob. 4. — The internal diameter of a steel cylinder is 20 cm. ; 
the wall is 10 cm. thick; the maximum safe permissible stress t 
in the material is 1200 kg. per sq. cm. Use this cylinder to make 
the foregoing theory clear to yourself. For example, (a) determine 
the permissible gas pressure when the first and the second approx- 
imations of the theory are used respectively; (6) assuming the 
cylinder to be just safe according to the first approximation, deter- 
mine the external diameter and the extra weight of steel, in per 
cent, necessary to make the wall sufficiently thick to be safe accord- 
ing to the second approximation; (c) make a drawing of the cross 
section of the cylinder to scale, and show graphically, by means of 
a curve, the distribution of stresses in the wall, according to both 

* When the student comes to the fourth approximation he will 
see that the effect of the lateral compression is not taken quite cor- 
rectly into account in Eq. (12a). However, for a thin sphere the result 
(13), based upon Eq. (12a), becomes identical with Eq. (5a). This 
fact shows the reasonableness of assumption (12), at least as an ap- 
I proximation. 
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approximations. Perfonn all computations very carefully, and 
make the sketch to a large scale, because the results are later to 
be supplemented with those obtained according to more accurate 
approximations. 

Prob. 5. — Plot a curve giving the ratio of the external to the 
internal diameter of forged steel guns for explosion pressures rang- 
ing from 3000 to 7000 kg. per sq. cm. The tensile stress in the 
material is allowed to go almost to the elastic limit, which in this 
case is approximately 3500 kg. per sq. cm. Draw the curve to a 
large scale because other curves are to be drawn according to more 
accurate approximations on the same sheet. 

Prob. 6. — Check the integration which leads to Eq. (13). 

Prob. 7. — r Solve Prob. 4 for a sphere, using the same data. 

Prob. 8. — Modify formula (10) for the case in which the ex- 
ternal pressure cannot be neglected. 

Prob. 9. — Modify formula (13) for the case when the external 
pressure cannot be neglected. 

Prob. 10. — Prove that for a very thin spherical shell, formula 
(13) becomes identical with Eq. (5a). 

Third Approximation. 

The Volume of the Wall is Assumed to Remain 

Unchanged under Stress. 

One obvious criticism of assumption (8) is that it 
indicates an increase in the volmne of the walls when 
a cylinder is subjected to an internal pressure. The 
original volume of the wall of a cylinder of length I is 

V ^tI (ri2 - ro2), 

while the volmne of the same cylinder imder pressure, 
according to the second approximation, is 

7' = tI [(ri + uy - (ro + u)^]. 

The difference 

V "V = Tl'2u(ri-ro) 

is positive, so that there would be an increase in volume 
if u were constant. It is much more reasonable, how- 
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ever, to believe that the total volume of the wall de- 
creases under pressure, the mner layers being pressed 
against the outer ones. The actual decrease in vol- 
ume not being known, the next best thing to assume 
is that there is no increase in volume. This is equiva- 
lent to saying that the tension t is distributed still less 
uniformly than in the second approximation, the values 
of u being larger at the inner layers and smaller at the 
external layers, so that the mside annuU press upon 
the outside ones. 

(a) Cylinder. — Let the values of u at the inside and 
the outside surfaces be tto and Ui respectively. Then 
the condition that there is no increase in volume becomes 

7' - 7 = ttZ [{n + ui^ - (ro + UoY] - tI (ri^ - ro^) = 0, 

or, after simpUfication, 

2 nui +Ui^ = 2 roih + U(?. 

The elongations Ui and Uq being very small compared 
with the radii of the cylinder, their squares may safely 
be neglected, and we finally get 

r\Ui = TqUq. 

This means that the extreme elongations are inversely 
as the radii of the cyhnder. It is reasonable to as- 
smne (though this is not logically necessary) that the 
same law holds inside the wall, or 

TiUi = xu = roVa = Const.* . . . (14) 

* Any function of the form 



I ^ ^ -^T^ + (x - fo) in - x)fix) 

I X 



satisfies condition (14) for x = n and x — r©, provided that / (x) does 
not become infinite at these values. We merely make the simplest 
assumption that/ {x) = 0. 
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Substituting the value of u from this expression into 
the general formula (7a) for t, we get 

t^E"^, .... (14a) 
and Eq. (1) gives 

Euol ) = po. 

Eliminating the product Euo from this formula by 
means of Eq. (14a) we finally get 

« = Po£»Y--^ (IS) 

\x/ ri — r© 

The maxunum tension is at the mner surface, and is 
equal to 

7*1 
fmax = P0- (16) 

n —To 

From this expression the necessary thickness of the 
wall for a given pressure and with a given inner diam- 
eter may be calculated. 

(6) Sphere. — The condition for no increase in vol- 
ume is 

r -.V=iT[(n+u,y- (ro +Uoy]-Uiri'-ro') = 0, 

or, neglecting again the higher powers of Ui and Uq, 

nhLi = roHtQ. 

Assuming, as in the case of the cyhnder, that the same 
law holds good within the wall, we obtain the relation 

riHii = xHl = ToHi^ = Const. . . . (17) 
Hence, Eq. (12a) becomes 

t = E^-^, (18) 
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and Eq. (2) becomes, after integration, 

^.2ro^(i-i)=|poro^. . . (19) 

Substituting the value of Eilq from Eq. (18) into Eq. 
(19) we finally get 

a formula similar to Eq. (15) for the cylinder. 

Prob. II. — Continue Prob. 4, using the third approximation; 
answer questions (a), (6), and (c). Supplement the drawing by 
showing the external diameter of the cylinder required, according 
to the three approximations, for the same gas pressure. Plot a 
curve of distribution of u knowing that for steel the modulus of 
elasticity E is equal to about 2 million kg. per sq. cm. 

Prob. 12. — In the first approximation t is independent of x, 
in the second approximation t is inversely proportional to x, in the 
third approximation t is inversely proportional to x*. Assuming 
that there are physical reasons why t should vary inversely as the 
nth power of x, n being a known number, write an expression for 
^maz similar to Eq. (16). 

Prob. 13. — For some materials (for instance, for cast iron), 
elongations Al/l in formula (6) are not proportional to tensions <, 
but increase more rapidly; in other words, the material is relatively 
weaker for greater stresses. Assume that the relation between the 
tension and the elongation may be represented by an empirical 
formula 



t 



^■(f)-''-(f)'. 



where E' and ^" are given constants, and deduce an expression for 
^max. Prove that with such a material, tensions are distributed more 
uniformly than with one that has a constant modulus of elasticity E. 

Prob. 14. — Modify formulas (15) and (20) for the case in which 
the external pressure cannot be neglected. 

Prob. 15. — Continue Prob. 7 for the sphere, according to the 
requirements stated in Prob. 11. 
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Prob. 1 6. — A set of spherical steel shells are to be cast for 
calorimetric tests of certain explosives. Show that these shells 
can be designed according to the formula 

J_ ^ Po 

Do 4:tmta — 2po' 

where 5 is the thickness of the shell, Do the inside diameter, and 
^maz is the allowable working stress in the material. 

Prob. 17. — Referring to Prob. 5, plot on the same sheet a 
similar curve according to the third approximation. 

Fourth Approximation. 

Lateral Deformations Taken into Account. 

(a) Cylinder. — An infinitesimal shell of radius x 
and thickness dx (Fig. 17) is subjected to three kinds 
of stresses, at right angles to each other: 

(1) Tangential tension t; 

(2) Radial compression p ; 

(3) Tension q parallel to the axis of the cylinder. 

The tangential tension resists the bursting of the 
cylinder and the radial compression is due to the inner 
layers of the wall pressing upon the outer layers under 
the influence of the gas contained in the cylinder. 
The longitudinal tension q is caused by the gas pressure 
on the enclosed ends of the cylinder; this is a stress 
which has heretofore been neglected. At a reasonable 
distance from the ends, stress q may be considered to 
be uniformly distributed over the cross section of the 
wall, so that 



PdTo^ - Piri^ 
ri2 - To^ 



g=gave- ^?"r;^ . . . . (21) 



Two more equations are necessary in order to de- 
termine t and p as fimctions of the radius x. One 
equation is obtained by considering the infinitesimal 
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shell of thickness da; as a thin cylinder subjected to 
pressure p on the inside, and to opposing pressure 
p + dp on the outside. Eq. (4) appUes to this case 
exactly, because t may be considered constant over the 
infinitesimal thickness dx. We thus have 

tdx = px — (p + dp) {x + dx), 

or, neglecting the infinitesimal quantity dp dx of the 

second order, 

tdx = -d{px) (22) 

A third relation between p, g, and t is obtained from 
the so-called law of lateral deformations. Namely, 
experiment shows that when an elastic body, say a 
rod, is subjected to a longitudinal tension t, it not only 
expands in the direction of t, but it also contracts 
in planes perpendicular to t This contraction, not 
caused directly by external forces, is called lateral 
contraction. Similarly, if the rod were compressed 
by an external force, the lateral deformation would be 
an expansion. Experiment further shows that lateral 
deformations are proportional to direct deformations, 
so that if a stress t produces an extension Ai/Z, the 
lateral contraction is equal to 



m\ir 



m 

where m is'a numerical coefficient. In this coimtry, 
this coefficient is usually assmned to be equal to 3. In 
France, where this coefficient is known as Poisson's 
ratio, it is taken equal to 4; C. Bach in Germany uses 
m = 10/3. Another way of stating the preceding law 
is that when a stress t is produced in a certain di- 
rection, stresses equal to —t/m are produced at right 
angles to t. 
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Thus, applying this law to our case, we can say that 
the longitudinal expansion of the cyUnder in the direc- 
tion of its axis is an algebraic result (a) of the direct 
tension produced by the pressures on its ends; (6) of 
the lateral contraction caused by tension t, and (c) 
of the lateral expansion produced by compression p. 
In other words, the annulus elongates in the direction 
of the axis of the cylinder as if an eflfective unit tension 
geff were applied to it, where 

m m 

In this equation we shall now make a simplifying 
assumption, namely, that the plane cross sections of 
the cyhnder remain planes after deformation. This 
condition is fulfilled in practice, since the ends of the 
cyUnder are sufficiently stiff to withstand any appre- 
ciable change in form. In other words, ?«« is inde- 
pendent of Xf being the same for all annuU. But, as 
has been stated above, the applied longitudinal stress, 
g, is also independent of x, at least at a reasonable dis- 
tance from the ends of the cylinder. Hence, accord- 
ing to Eq. (23), {t — p) must also be independent of 
X, or 

^ - p = 2 A = Const., . . . (24) 

where 2 A is a constant to be determined later. 

From Eqs. (21), (22), and (24) the three unknown 
stresses p, g, and t may be calculated. Of these, q is 
directly known from Eq. (21); so that Eqs. (22) and 
(24) may be solved together for p and t. Substituting 
the value of t from Eq. (24) into Eq. (22) gives 

(p+2A)dx=^ —(pdx + xdp)f 
or 

'-2Adx = 2pdx+xdp. 
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This differential equation is easily integrated by multi- 
plying both sides by x; it is reduced then to 

—2 Ax dx = d {x^)y 

or, after integration, 

-Ax'^ +B = x^, 

where B is the constant of integration. Solving for 
p gives 

P=|-^ (25) 

Hence, from Eq. (24), 

t = p + 2A ==~ + A. . . . (26) 

x^ 

The constants A and B are determined from.Eq. (25), 
knowing that p = Po inside the cyUnder and p = pi 
outside the cylinder, so that 



B A 

B A 



(27) 



Solving these equations for A and B, it will be found 
that both A and B are positive. Hence, a comparison 
of Eqs. (25) and (26) will easily show that the tan- 
gential tension t is always munerically greater than the 
radial compression p. Both reach their maximum for 
X = ro, viz., at the inner surface of the cyUnder. 

(6) Sphere. — Consider again the infinitesimal shell 
of diameter x and of thickness dx. The shell is in 
equiUbrimn under the influence of pressure p on its 
inner surface, pressure (p + dp) on the outside surface, 
and tension L Applying Eq. (5) to this case, we get 

2txdx = px^ — (p + dp) {x + dx)% 






. - '^t 
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or, performing the multipUcation and neglecting the 
infinitesimal terms of higher order, 

2tedx = -d{px^) (28) 

This equation takes the place of Eq. (22) for the cylin- 
der. The law of lateral deformations fmnishes two 
equations for the sphere: 

^,ff = ^H =<«! + £; . . . (29) 

X mm 

To deduce Eq. (29) consider a small square portion of 
the infinitesimal shell of thickness dx. This portion 
expands in two perpendicular directions under the m- 
fluence of the tensions f, and is compressed radially by 
the pressure p. Thus, according to the law of lateral 
deformations, 

m m 
But, according to the fundamental Eq. (7a), 

^eflf = E — f 
X 

from which Eq. (29) follows. By comparing Eq. (29) 
with Eq. (12a) it will be seen that if p could be neglected 
and if m were equal to 2, the two equations would be 
identical. This shows the error, or rather the assmnp- 
tions, made in Eq. (12a). It is proper to use Eq. (7a) 
for the sphere, because it is applied here to the effective 
tension, the lateral deformations having been properly 
accoimted for. Eq. (30) expresses the fact that the 
. radial pressiu'e p is further increased by the lateral 
:y:f J)ressures t/m due to the tangential tensions t in two 



V. 
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mutually perpendicular directions. Thus the effective 
radial pressure is 

Peflf = P H 1 

mm 

On the other hand, it is proved below that du/dx 
represents the relative radial compression which in the 
case of dilatation is denoted by Al/l in Eq. (6). Hence, 

and Eq. (30) follows immediately. 

To prove that du/dx represents the relative contrac- 
tion, consider a point a on the inside surface of the 
infinitesimal shell, and a point &, on the same radius, 
but on the outside surface of the same shell. The 
origmal distance between the points is dx. After the 
deformation, the distance of a from the center be- 
comes X +u, while that of & becomes (x + dx) + 
(u + du). Now 

Ai new length ab — original length ab 
I " original length a6 

or 

AZ __ { [(x + dx) +(u + du)] — (x + u)\'-dx _ du 
I dx dx 

which was to be proved. The minus sign in Eq. (30) 
is necessary before du/dXy because u decreases from 
the center to the periphery and consequently du/dx is 
negative, while the right-hand side of the equation is 
essentially positive. 

The three equations (28), (29), and (30) contain 
three unknown functions of x, namely p, t, and u. To 
solve these equations we substitute the values of t 
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and p from Eqs, (29) and (30) into Eq. (28). This 
gives a differential equation between u and x, so that 
u can be determined as a function of x. Substituting 
this value of ti in Eqs. (29) and (30) p and t are ex- 
pressed as functions of x, and the problem is thus 
solved. 

Solving (29) and (30) for p and t gives 

Substituting these values in Eq. (28) we obtain^ after 
reduction, 

x^p{ + 2x^^2u (33) 

dx^ dx ^ ^ 

This equation may be integrated according to the gen- 
eral rules for the integration of differential equations 
of this type,* but it can also be solved as follows: 
Multiply both sides by x and add to both sides x^ du/dx; 
the result is 

dx»^ dx ^ dx 



or 



dx\ dxl dz 



Integrating both sides, we get 



a^^ = x*u-3C„ .... (34) 



* See Probs. 26 and 26 below. 
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where 3 C2 is a constant of integration. To integrate 

again, we divide both sides of the equation by x*; the 

result is 

du 

dx 3 C2 



or 



x^ X* 



\X/ X* ' 



da; 
Integrating both sides, we find 

u ^ C2 

x ar 

where Ci is another constant of integration. Solving 
for u gives 



x^ 



u = CiX+^^ (35) 



This is the expression f or 2i as a function of x. Substi- 
tuting it in Eqs. (31) and (32) we find 

< = ^S-^-& + -^-^ii' • • (36) 
(m + 1 XT m — 2 ) 



P = ^j 



m + 1 x^ m — 2 



The constants C2 and Ci are determined by applying 
Eq. (37) to the terminal conditions, viz., p = po when 
X = ro, and p = pi when x = n. 

Pxob. 18. — Continue Probs. 4 and 11, using the formulas of the 
fourth approximation. 

Prob. 19. — Prove that A =^ q,ao that qeB is never greater than 
from Jtojg, ifwi = 3to4. 

Prob. 20. ^— Continue Probs. 5 and 17, using the formulae of 
the fourth approximation. 
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Prob. ai. — From Eqs. (36) and (37) it follows that 2 « — p = 
Const. This corresponds to Eq. (24) for the cylinder. Assuming 
this relation to be given, what is the simplest way to deduce Eqs. 
(36) and (37)? 

Prob. 22. — Represent Eqs. (36) and (37) in their final practical 
form, without Ci and Cj. 

Prob. 23. — Show that the formula given in Prob. 16 should be 
modified to 



— = 1 /\/ ^ax + Po _ i\ 
Do 2\^tmBx-po / 



Prob. 24. — Continue Probs. 7 and 15, using the formulse of the 
fourth approximation. 

Prob. 25. — Solve Eq. (33) by first reducing it to one with con- 
43tant coefficients, and then solving the latter by the usual method 
of the auxiliary equation. Hint. — Use the substitution x = e'. 

Prob. 26. — Solve Eq. (33) by assuming a priori that the solu- 
tion is of the form 



Fifth Approximatioii. 
True Deformation as a Mea^sure of Effective Stress. 

(a) Cylinder. — 'Eqs. (21), (25), and (26) give the 
values of the component stresses at any point of the 
cylinder wall. Under the influence of these stresses 
the material is strained, each component stress pro- 
ducing a direct deformation in its own direction and 
lateral deformations of the opposite sign in perpendicu- 
lar planes. Therefore, the resultant actual deforma- 
tions are not proportional to the component stresses, 
but must be calculated from the component stresses 
according to the law of lateral deformations. It is not 
the stress, but the actual deformation that determines 
the safety of the material at a given point. 

Let, for instance, a cube whose edge is one cm., be 
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subjected to a compression of 100 kg. upon two oppo- 
site sides, and to a tension of 60 kg. upon two other 
opposite sides; let the third pair of sides have no 
forces appUed to them. The total deformation in the 
direction normal to the first pair of sides corresponds 
not to a stress of 100 kg. per sq. cm., but to a stress 
of 100 + J X 60 = 115 kg. compression per sq. cm., 
assuming m = 4. Similarly, the deformation normal 
to the second pair of sides corresponds to a stress of 
60 + J X 100 = 85 kg. per sq. cm., while in the direc- 
tion perpendicular to the third pair of sides on which 
no external forces are applied, the deformation corre- 
sponds to a tensile stress of J X 100 — J X 60 = 10 kg. 
per sq. cm. 

Thus, while expressions (21), (25), and (26) are cor- 
rect in themselves, the values of p, g, and t calculated 
from these equations do not give the true stresses in 
the material. The true or the effective stress parallel 
to the axis of the cylinder is not q, but is geff as given 
by Eq. (23). Similarly, the true tangential stress is 
not f, but 

te« = t-^+^, .... (38) 
m m 

and the true radial compression is 

m m 

The effective values may be smaller or greater than the 
component stresses p and t, according to the values of 
the radii of the cylinder and of the applied pressiu'es. 
Substituting the values of t and p from Eqs. (25) and 
(26) mto Eqs. (23), (38), and (39), and remembering 
that A ^ q (see Prob. 19), we get 



78 



BNGINEERINO MATHEMATICS. (Chap. VL 



ffrf-? 



p«ir = 



4« = 










• 



. (40) 



The value of <«« at the inner surface of the cylinder 
usually determines the required thickness for a given 
factor of safety, or vice versa. 

(6) Sphere. — Instead of being given by Eq. (36), 
the true tangential stress, according to Eqs. (29) and 
(35), is equal to 

«.. = ^(| + Ci) (41) 



Similarly 



Pes 



=^(^'"4 • • • (^) 



The constants Ci and C2 are determined, as before, 
by applying Eq. (37) to the terminal conditions. 

Prob. 27. — Continue Probs. 4, 11, and 18, using the formulfie 
of the fifth approximation. 

Prob. 28. — Write down fonnulse for the true stresses in the 
material of a cylinder with very thin walls. 

Prob. 29. — Deduce formulae for the strength of a thick cylinder 
subjected to an external pressure, for example, a fire tube in a 
steam boiler, or a pipe culvert. 

Prob. 30. — A compoimd gun is to be formed out of two con- 
centric cylinders in the following way: The inside radius of the 
inner cylinder is to be ro, the outside radius p. The inside radius 
of the outer cylinder is left less than p by a small amount a; the 
outside radius is to be n. The outer cylinder (sometimes called 
the jacket or the hoop) is heated so that its inside radius becomes 
equal to p, and then it is shrunk on the inner cylinder. Assuming 
fo, p, and n to be known, show how to determine a so that during 
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firing when a pressure po is developed inside, the effective tensile 
stresses on the inner surfaces of both cylinders will reach the elastic 
limit of the material. 

Note. — The purpose of using hoops on guns is to make them 
lighter. In an ordinary gun the material is worked up to the safe 
limit only on the inner surface of the cylinder. The hoop gives 
an original compression to the inner cylinder, so that the explosion 
of the gases must first overcome this compression before it pro- 
duces a tensile stress t. Again, the inner surface of the hoop may 
be worked up to the elastic limit, so that the material is utilized 
better than in a solid gun. 

Piob. 31. — Determine the neces&ary thickness of the barrel of 
a gun whose internal diameter is 15 cm., and in which the pressure 
during firing reaches 4000 kg. per sq. cm. Then design the same 
gun as a compound cylinder with a hoop, selecting the radii of the 
internal cylinder and of the hoop so as to use a minimum of mate- 
rial. See Prob. 5 for the permissible stress. 

Prob. 32. — Continue Probs. 5, 17, and 20, using the formula 
of the fifth approximation. 

Prob. 33. — Continue Probs. 7, 15, and 24, using the formute 
of the fifth approximation. 

Prob. 34. — Write down formulae for true stresses and for def- 
ormation u for a very thin elastic sphere. 

Prob. 35. — Write down formulsB (41) and (42) in their final 
practical form, with the constants Ci and C2 expressed through 
the known quantities. 

Prob. 36. — The above given theory of thick cylinders presup- 
poses that longitudinal tensions q are independent of x. This 
does not hold true near the ends of the cylinder. If a cylinder has, 
say, hemispherical ends, the stresses in the connecting plane come 
out different according to whether they are calculated from the 
cylindrical or the hemispherical side. It would be a suitable prob- 
lem for original research to develop formuke for stresses in that 
part where the cylinder connects with the hemisphere, so that the 
changes in the values of the stresses would come out to be gradual, 
as they are in reality, and not sudden. 
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